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We report on the feasibility of measuring the top Yukawa coupling in the process: 
e + e~ —¥ ttH. This measurement is crucial to test the mass generation mechanism for 
matter particles. Since the cross section for this process attains its maximum around 
\fs — 700 GeV, most of the past studies were done assuming this energy region. It 
has been pointed out, however, that the QCD threshold correction enhances the cross 
section significantly and might enable its measurement at y/s = 500 GeV, which will 
be accessible already in the first phase of the ILC project. We have implemented 
this threshold enhancement into our ttH event generator and carried out Monte Carlo 
simulations. Our results show that ttH events can be observed with a significance 
of 4.1cr with no beam polarization and 5.4 a with the e~ and e + beam polarization 
combination: (—0.8, +0.3). 



1 Introduction 

The standard model of elementary particle physics is based on two pillars: one is the gauge 
principle and the other is the electroweak symmetry breaking and mass generation mech- 
anism. The first pillar, the gauge principle, has been tested by precision electroweak mea- 
surements. On the other hand, the second pillar has not yet been tested. In order to confirm 
this second pillar we have to measure the Higgs self-coupling and the top Yukawa coupling. 

In this study, we investigate the feasibility of measuring the top Yukawa coupling at 
500 GeV with the process: e + e~ — > ttH. Since the top quark is the heaviest among all 
the matter particles, the measurement of its Yukawa coupling will be the most decisive 
test of the mass generation mechanism for matter particles. Since the cross section for the 
e + e~ — > ttH process is 2-3 fb even near its maximum reached at around \fs = 700 GeV, 
most of the past studies assumed the measurement energy in this region (|T],j2]). It has 
been pointed out, however, that the QCD threshold correction enhances the cross section 
significantly and might open up the possibility of measuring the top Yukawa coupling at 
sfs = 500 GeV, which is within the scope of the first phase of the ILC project ([3] -[19]). In 
order to investigate this possibility we have implemented this threshold enhancement into 
our ttH event generator and carried out Monte Carlo simulations. 

In the next section we begin with clarifying the signatures of the ttH production and 
list up possible background processes that might mimic the signal. We then describe our 
analysis framework used for event generations and detector simulations in section 3. The 



event selection procedure for the generated events is elaborated in section 4, considering 
characteristic features of the background processes. The results of the event selection are 
given in section 5. Section 6 summarizes our results and concludes this report. 

2 Signal and Possible Background 

The Feynman diagrams for the e + e~ — > tiH process followed by t(i) — > b(b)W decays are 
shown in Figure [U Notice that the first and second diagrams contain the top Yukawa 
coupling, which we want to measure. The signatures of tiH events depend on how the 
H and the Ws decay. In this study we concentrate on the dominant decay mode: H — > 
bb(Q8%). The signal events hence have four b jets and two Ws. The ttH events can then 
be classified into 3 groups (8-jet, l-lcpton+6-jet, and 2-lepton+4-jet modes) corresponding 
to the combinations of lcptonic and hadronic decays of the two Ws. For W^s that decayed 
lcptonically we cannot reconstruct their invariant masses due to missing neutrinos. On the 
other hand, for the Ws that decayed hadronically we can reconstruct their masses and use 
them as a signature. For the t or the i with a hadronically decayed W we can also use the 
invariant mass of the 3-jet system to test if it is consistent with the top mass. 




Figure 1: Fcynmann diagrams for the ttH process 

Possible background processes that might mimic the signatures of the ttH production 
include e + e~ — > tiZ , ti, and tig followed by g — > bb. The cross sections for these background 
processes are plotted in Figj2] together with that of the signal. Notice the smallness of 
the contribution from the third diagram in Figfl] which does not contain the top Yukawa 
coupling. We can hence determine the top Yukawa coupling by just counting the number 
of signal events unless they are swamped by the background; the signal cross section is only 
~ 0.5 fb with no beam polarization. 

The production cross section for the tiZ background is 1.3 with no beam polarization. 
It has four 6-jets and two Ws in the final state just like the signal, if the Z boson decays into 
bb (15%). In this case the only difference that one can tell on an event-by-event basis lies 
in the invariant mass of the bb system, which should be consistent with M# for the signal 
and Mz for the background. The ti production, on the other hand, has only two 6-jets 
in the final state. If reconstructed correctly, it could not be the background. Since the ti 
production cross section (~ 500 fb) is much larger than that of the signal, however, a small 
fraction of mis-reconstruction or failure in 6-tagging may lead to significant background 
contamination. The tig production followed by g — > bb decay has the same signatures as the 



* This value is with QCD threshold enhancement similar to that expected for the signal process. Without 
the correction the cross section is 0.7 fb. 



with no beam polarization 
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Figure 2: Production cross section of the signal, tiH, together with those of the main 
background processes, ttH,ttZ,ti,ttg, as a function of the center of mass energy for no 
beam polarization. 



signal in terms of the number of 6-jets and the number of Ws. As with the ttZ background 
the only difference is the invariant mass of the bb system. Its production cross section is also 
of the same order, 0.7 fb, as that of the ttZ background. 



3 Analysis Framework 



For Monte Carlo simulations, we generated signal and background events by using an event 
generator package (physsim[20]), which is based on full helicity amplitudes calculated with 
HELASjH] including gauge boson decays, thereby correctly taking into account angular 
distributions of the decay products. The 4-momcnta of the final-state quarks and leptons 
were passed to Pythia6.4[22 for parton showering and hadronization. The resultant particles 
were then swum through a detector model (see Table [T] for detector parameters) defined 
in our fast Monte Carlo detector simulator (QuickSim[53]). In the event generations we 
used a(M z ) = 1/128, sin 2 9 W = 0.230, a s = 0.120, M w = 80.0 GeV, M z = 91.18 GeV, 
Mt = 175 GeV, and Mr = 120 GeV. We have included the initial state radiation and 
beamstrahlung in the event generations. The unique point of this study is the inclusion of 
the QCD threshold enhancement to the ti system (see Figl3]) for the signal event generation, 
which plays an important role especially in a low energy experiment: about a factor of 2 
enhancement at ^/s = 500 GeV. 
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Figure 3: Invariant mass distribution for the tt sub-system. 



Table 1: Detector Parameters, where p,pr and E are measured in units of GeV 



Detector 


Performance 


Coverage 


Vertex detector 


cr fc = 7.0 8 (20.0/p) / sin 3/2 9fim 


|cos0| < 0.90 


Central drift chamber 


<j Pt /P t = 1.1 x IQ-^pt 0.1% 


| cos6>| < 0.95 


EM calorimeter 


a E /E = 15%/VE®l% 


|cos0| < 0.90 


Hadron calorimeter 


ct e /E = 40%/Ve®2% 


|cos0| < 0.90 



4 Event Selection 

4. a Definition of our signal (l-lepton+6-jet mode on tiH) 

As explained in section [2] we can classify the tiH signal events into the following three decay 
modes according to how the two Ws from t and i decay: 

1. 8-jet mode (45%) 

2. 1-lcpton + 6-jet mode (35%) 

3. 2-lepton + 4-jet mode (7%) 

where the lepton is required to be either e^- or ^ and the final-state H to decay into the 
dominant bb state. Notice that in all of these three modes we have four 6-jets in the final 
states, which makes the separation of the ti background easier. In this study we concentrate 
on the 1-lepton + 6-jct mode as our first step because the branching ratio is not so low and 
the number of jets is not so high. 

As shown in FigsfJ]and[S]tlic signatures of our signal are 
• an isolated energetic e^ 1 or /z^, 



• six jets including four 6-jets, two of which form a H boson, 

• the remaining two jets being consistent with a W boson, and 

• one of the two unused 6-jets together with this W candidate comprising a t quark. 




b isolated charged 
lepton 



Figure 4: Schematic diagram 
defining our signal signatures 



Figure 5: Invariant mass distributions for the hadronically 
decayed W, t, and H, which are reconstructed using gen- 
erator information. 



In what follows we will elaborate selection cuts designed to single out these signatures. 



4-a.l Isolated lepton search 

Our event selection starts with the search for a lepton coming from a W — > Iv decay. Such 
a lepton from W tends to be energetic and isolated from the other tracks. In order to find 
such an isolated lepton, we consider a cone around each lepton track (see Fig|6|) and define 
the cone energy to be the sum of the energies of the other tracks in the cone. Figure [7] 
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Figure 6: A cone around lepton track Fi S ure 7: Cone ener §y dis tribution o f isolated lep- 
ton: cut boundary y = y // 6(x — 15) 

plots the cone energy against the lepton energy. The energetic isolated leptons from Ws 
have to have a high lepton energy and a low cone energy, hence populating the bottom edge 
region (black points), while leptons from heavy flavor jets are likely to be less energetic and 
have a higher cone energy (gray points). The smooth curve in the figure is our cut to select 
energetic isolated leptons. 



4-.a.2 Forced 6-Jet clustering 



After finding and eliminating an energetic isolated lepton, we perform jet clustering to make 
six jets. For the jet clustering we use a variable Y defined by 

y M U 

E 2 

visible 

We keep putting tracks together to form a jet while Y < Y cu t- By adjusting the Y cu .t value, 
we can make arbitrary number of jets. We hence force the events to cluster into six jets by 
choosing an appropriate Y cu t value on the event- by- event basis (forced 6-jet clustering). 

4. a. 3 Y cut cut 

The Y cut value for a ti background event to form six jets should be lower than the one for 
a signal tiH — > tibb event because, after the energetic isolated lepton requirement, the ti 
event can hardly have more than four jets. The difference in the Y cut value distributions 
between UH(H — > bb) and ti is shown in FigJHl As seen in the figure, by cutting Y cut values 
at 0.002 we can reduce the ti background effectively. 
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Figure 8: Ycut value distribution after isolated lepton finding 



4-a.4 mass cut 

After performing the jet clustering, we try to identify which jet is coming from which parent 
parton. We want to separate the correct combination from the other combinatorial back- 
ground. Mass cut comes in handy to reduce the combinatorial background. Looping over 
all the 2-jet combinations we look for a pair having an invariant mass within the window 
of ±15GeV from the nominal W mass of 80.0 GeV. From the remaining four jets we pick 
up one and attach it to the just found pair making a W candidate to see if the resultant 
3-jet system has an invariant mass within ±25 GeV from the nominal t mass of 175 GeV. 
If it does we search for a pair from the three jets left over that is within the mass window 
of ±15 GeV from the nominal H mass of 120 GeV. Since these mass cuts are rather loose 
there is a significant chance to have multiple combinations that pass them. For such a case 
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we define a x 2 variable with 



' A/ 2-jct(w) ~ M w ' 

0~Mw 



M 3-jet(t/t) ~ M * ' 
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and select the combination with the smallest x 2 value. Fig|9] shows the mass distributions 
for the best combinations. Although W and t/i peaks are present for both the signal and 
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Figure 9: Invariant mass distributions after the cut on Y cu t values. Black open histograms 
are for the signal and gray histograms are for the ti background. 



the tt background, a H peak is seen only for the signal process. The H peak is, however, 
swamped in the tt background. 

4- a. 5 b-tagging by the n-sig. method 

For the tt background rejection, 6-tagging is very powerful since the signal ttH(H — > bb) 
process has four 6-jets, while the tt background process has only two 6-jets. For 6-tagging 
we use the so called n-sig. method descrived as follows. 

Figure [TOl sketches a jet from the interaction point (IP), which includes a &-hadron. The 
6-hadron decays at distance from the IP due to its long-life. It makes the 6-jet to have some 
tracks which are away from the IP. When the distance (£) between the IP and a track is 
larger than a given value {mat, ), the track is defined as an off- vertex track. A jet is recognized 
as a 6-jet if the number of such significantly off- vertex tracks exceeds a certain cut value 
(n). In this analysis, we define tight 6-tagging with a tagging condition: (m, n) = (3.0,2) 
and loose 6-tagging with (m, n) = (2.0,2), and require all of the four &-jet candidates have 
to satisfy the loose ^-tagging condition and there has to be at least one tight fr-tagged jet 
from each of the H and t/i candidates. 



Figure 10: n-sig. method 



The mass distributions after the 6-tagging are shown in FigfTTJ We can see that the tt 
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Figure 11: Invariant mass distribution after using both Y cut and b-tagging 



background has been suppressed effectively. As mentioned above the ttZ and tig (g — > bb) 
background events have similar signatures as a signal and can be separated only with the 
invariant mass of the H candidate. In the next section we summarize the results of our 
event selection including these remaining background processes. 

5 Results 

In order to estimate the feasibility of measuring the top Yukawa coupling we need to specify 
the beam polarization and the integrated luminosity. In this study we assume an integrated 
luminosity of 1 ab _1 . As for the beam polarization, it is worth noting that only the left-right 
or right-left combination contributes to the signal and background cross sections because of 
the 7 M coupling of the beam particles to the vector bosons (f/Z) in the intermediate states. 
It is hence sufficient to know the cross sections for the beam polarization combinations: 
(e~,e + ) = (— 1, +1), (+1, — 1). Tabl(J2] shows these cross sections. 

For both of the beam polarization combinations: (— 1,+1) and (+1, —1), we have gen- 
erated 50k events each for the ttH, tiZ, and tig (g — > bb) processes, and 5M events for 
the tt background. We performed the event selection described in the previous section and 
tabulated the results normalized to an integrated luminosity 1 ab _1 in Table[3] assuming the 
cross section shown in Tabled 

The corresponding distributions for the reconstructed W, t/i, and H candidates are 
shown in FigfT2]for the beam polarization combination: (—0.8, +0.3). We can see a clear 



Table 2: Cross sections at y/s = 500 GeV. ttH and ttZ are with QCD threshold enhance- 
ment. (-1, + !)/(+!, -1) corresponds to (e£,e^)/(e^,ej), respectively. 



Beam Polarization 


(-V 


H) 


(+1, 


■1) 


ttH 


1.24 


[fb] 


0.540 


[fb] 


ttz 


2.18 


[fb] 


0.712 


[fb] 


tt 


720. 


[fb] 


309. 


[fb] 


tig (g -> 66) 


1.93 


[fb] 


0.859 


[fb] 



Table 3: Cut Statistics (normalized to 1 ab 1 ) 



Beam Polarization 




(0.0,0.0) 






(-0.8,+0.3) 




Processes 


ttH 


ttz 


tt 


tig (66) 


ttH 


ttz 


tt 


ttg (66) 


No Cut 


449.0 


1340.0 


514040.5 


697.5 


759.0 


2407 


863500.4 


1159.6 


so. 1 ep 1 


159.4 


435.9 


209718.4 


242.2 


269.4 


783.0 


303879.0 


397.7 


Ycut (6 jets) > 0.002 


139.2 


307.8 


22851.3 


152.5 


235.4 


552.9 


38477.2 


249.6 


btag & mass cut 


23.0 


12.2 


11.9 


6.9 


38.9 


21.8 


19.7 


11.3 



evidence of signal events over the background in each of the three mass distributions. 




Figure 12: Mass distributions (cumulative) for the final selected sample for the beam 
polarization combination: (—0.8, +0.3). 



In the case of no beam polarization 23.0 signal events are left with 31.0 background events 
total. On the other hand we have 38.9 signal events with 52.8 background events total at the 
end of the event selection. The signal significance is 4.1 a for the polarization combination: 
(0,0) and 5.4cr for the polarization combination: (—0.8, +0.3). Since the number of the 
signal events is proportional to the square of the top Yukawa coupling (gy), we can easily 
translate these numbers to its expected precisions: Agy/gy = ±0.12 and ±0.093 for the 
beam polarization polarization combinations: (0,0) and (—0.8, +0.3), respectively. 



6 Summary and Conclusion 



We have performed a feasibility study of measuring the top Yukawa coupling at s/s = 
500 GcV, taking advantage of the QCD threshold enhancement to the ti sub-system. For this 
study we have implemented the threshold enhancement in the tiH and tiZ event generators 
in the physsim package. It is found that for an integrated luminosity of lab -1 we can 
observe the tiH process with a significance of 4.1 a without beam polarization, and 5.4 a 
with the beam polarization combination: (e~,e + ) = (—0.8,-1-0.3). These numbers show that 
we can measure the top Yukawa coupling to an accuracy of about 10% at y/s = 500 GeV, 
which is the energy already available in the first stage of the ILC. 
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